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We present a class of cancellation mechanisms to suppress the total contributions of Barr-Zee 
diagrams to the electron electric dipole moment (eEDM). This class of mechanisms are of particular 
significance after the new eEDM upper limit, which strongly constrains the allowed magnitude of 
CP-violation in Higgs couplings and hence the feasibility of electroweak baryogenesis (EWBG), 
were released by the ACME collaboration in 2013. We point out: if both the CP-odd Higgs-photon- 
photon (Z boson) and the CP-odd Higgs-electron-positron couplings are turned on, a cancellation 
may occur either between the contributions of a CP-mixing Higgs boson, with the other Higgs 
bosons being decoupled, or between the contributions of a CP-even and a CP-odd Higgs bosons. 
With the assistance of the cancellation mechanisms, a large CP-phase in Higgs couplings with 
viable electroweak baryogenesis (EWBG) is still allowed. The reopened parameter regions would 
be probed by the future neutron, mercury EDM measurements, and direct measurements of Higgs 


CP-properties at the LHC and future colliders. 


INTRODUCTION 


The baryon asymmetry in the Universe (BAU) nowa- 
days, i.e., 

a ~ (0.7 — 0.9) x 1071 £ 0 (1) 
has puzzled people for more than half a century. Here 
s is entropy density of the Universe. Among various 
dynamical mechanisms to solve this puzzle, electroweak 
baryogenesis (EWBG) falls in the most popular class, due 
to its potential testability at the Large Hadron Collider 
and in the other experiments. A generic feature of the 
EWBG is that the CP phases employed to generate the 
cosmic baryon asymmetry need to enter the couplings 
between the Higgs sector and particles which either ex- 
ist in the Standard Model (SM) or are introduced in new 
physics, no matter the CP-phases are flavor-diagonal, off- 
diagonal [3], or flavor-decoupled. Otherwise, these CP- 
phases are decoupled from electroweak phase transition 
(EWPT) and the EWBG will never be implemented. The 
measurement of the Higgs CP-properties therefore pro- 
vides important information to solve the BAU puzzle. 

Motivated by this, the CP-properties of the Higgs bo- 
son discovered in 2012 [4] have been extensively stud- 
ied by both theorists [5}{9] and experimental groups [IJ] 
since its discovery [4], by using a method of direct mea- 
surements at the LHC. Given the limited statistics, how- 
ever, the sensitivity of the LHC at this stage is still low. 
On the other hand, fast progress has been made in in- 
direct measurements. Using the polar molecule thorium 
monoxide (ThO), the ACME collaboration reported an 
upper limit on the eEDM recently [12] [56] 


Ide] < 8.7 x 107-7%ecm (2) 


at 90% confidence level, an order of magnitude stronger 
than the previous best limit. This limit severely con- 
strains the allowed magnitude of CP-phases in the Higgs 


couplings [6H9] via Barr-Zee diagrams, causing a tension 
between the observation and the CP-phase required for 
successfully implementing EWBG (e.g., see [10] where 
the expected projection of the eEDM bounds to the 
EWBG in the MSSM was studied.) [57]. 

In this letter we point out that in these studies a cru- 
cial effect was more or less ignored, which can dramati- 
cally change the conclusions. This is due to the fact that 
generally both the CP-odd Higgs-photon-photon and the 
CP-odd Higgs-electron-positron couplings can be or tend 
to be turned on. These two couplings contribute to the 
eEDM separately and simultaneously. If there exists a 
cancellation between their contributions (as we will show 
below in two different contexts: the type-II two Higgs 
Doublet Model (2HDM) where the tree-level CP-phase 
arises from the pure Higgs sector, and the Minimal Su- 
persymmetric Standard Model (MSSM) where the tree- 
level CP-phase arises from Higgs-superparticle interac- 
tion sectors.), even if the magnitudes of the CP-phases 
in Higgs couplings are large, the current ACME bound 
can be well-satisfied. In such a case, EWBG can still be 
successfully implemented. 


GENERAL ANALYSIS 


In an effective Lagrangian for a Higgs sector, the rele- 
vant operators are given by: 
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where F y is field strength of photon, with Pu 
(1/2)euvpo F”, Vav is field strength of photon and 


boson, with Viv = (1/2)euvpo V®?, and 0 = tanm! 
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defines the CP-phase of the Yukawa couplings. These op- 
erators can be inserted in the Barr-Zee EDM diagrams. 
Integrating out the internal degrees of freedoms, we have 


Lef = —ideēo"” y5e0, Ay 5 (4) 
with its contribution to the eEDM given by 
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Here v =246 GeV is the normalized vacuum expectation 
value (VEV) of the Higgs fields, and A$;y (Äiv) is the 
relevant scale for the hi F”” Vay (hiF”” V u) operator. It 
is clear that the Barr-Zee contributions depend on not 
only the CP-odd Higgs di-photon coupling č}, but also 
the CP-even one c} if the Higgs bosons have a CP-odd 
coupling with electrons (č, Æ 0). 
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FIG. 1: Two cancellation mechanisms of the Bar-Zee contri- 
butions to the eEDM. Upper: cancellation occurs between the 
contributions of a CP-mixing Higgs boson. Bottom: cancel- 
lation occurs between the contributions of a CP-even and a 
CP-odd Higgs bosons. 


The ACME measurement greatly improved the current 
bound on the eEDM, leading to 


a Ny \ 86 Iog { Atv 
-c č l ech l .14 
X ceČ log ma + čec, log ma < 0.14 (6) 
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This strongly constrains the allowed CP-violation in a 
single Higgs coupling, e.g., in the case with one (SM-like) 
Higgs only and meanwhile če = 0 [13], unless this Higgs 
boson is decoupled. However, if a cancellation occurs 
among these inference terms, CP symmetry is allowed 
to be significantly violated, without contradicting with 
the current eEDM bound. Below we will use the type 


II 2HDM and the MSSM to show two different cancel- 
lation mechanisms, both of which are mainly motivated 
by EWBG: (1) cancellation occurs between the contribu- 
tions of a CP-mixing Higgs boson, while the other Higgs 
bosons are decoupled (see the upper diagrams in Fig. D; 
and (2) cancellation occurs between the contributions of 
a CP-even and a CP-odd Higgs bosons (see the bottom 
diagrams in Fig. (i). 


TYPE IT 2HDM 


As an illustration, we consider type II 2HDM with a 
softly Z2 symmetry (¢1 + —¢ı and ¢2 > ¢2) [14]. Its 
tree-level Higgs potential is given by: 


V = = (lor)? + 2 (hon)? + As(dl on) (her) (7) 
+ ra($4¢2)(¢h¢1) + : [stolga + h.c. 


— 5 {m2, (6dr) + [mfa(blo2) + b.c.] + maloto) 
Here m 2 and Às are complex parameters. Their relative 
phase Arg(A5m43) leads to CP violation in the Higgs sec- 
tor. We take the convention that both Higgs doublets 
¢1,2 carry a hypercharge of one unit and that the general 
Higgs VEVs are 


whe) @=(,e),  @ 


with sin? 8 = |vo|?/(|v1|? + |vel?), vı = v cos 8/V2 and 
|va| = vsin 8/V2. Then the unitary matrix R, defined to 
diagonalize the Higgs mass matrix M 


RMR” = diag( Mh, , Mi, Mhs)» (9) 


in the mass engenstate (h1, h2, h3) can be easily figured 
out, given by 
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with c; = cosa;, si = sina;. Here a, ay, Qe are mixing 
angles between two CP-even Higgs, the light CP-even 
and the CP-odd Higgs, and the heavy CP-even and the 
CP-odd Higgs, respectively. The angular range, beyond 
which R is repeated, can be chosen as 0 < a < T, =T < 
ay < T and —17/2 <a. < 7/2. 

The tree-level hı couplings rescaled by the SM values 
are given by 


COS Q COS Œp sin @ COS Ap 
wor sing ” ee cosB ’ 
Č = — cot sin Qp, Č = Če = — tan ĝ sin ap , 
ay = cosapsin(6 — a) (11) 


here hy is SM-like and ay represents the hıWW and 
h,ZZ couplings. The CP-phase of the top Yukawa cou- 
pling 6; is given by 
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c 
tan 6; = E tan ap . (12) 
These tree-level effective coupling further contribute to 


cy and é, at loop level 


a = QF cz /2 = 24/9 
č, = -3Q76,/4 = —č/3 
W = —Tay/8. (13) 


Note that the signs of c$, or ¢, are derived in the conven- 
tion of Eq. (3). To achieve the cancellation indicated by 
Eq. f5}, we have hW Zef Č ce ~tanB~ 1. 

In this setup, the lightest Higgs boson hı leads to a 
leading-order contribution to the eEDM ; 
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via the Barr-Zee diagrams [I8], where z = m?/mj,, 


zw = my/m;,, and the loop functions f(z) and g(z) 
are given in [18]. Numerically, we have f(z:) = 1.0 and 
g(zz) = 1.4. These quantities depend on three free pa- 
rameters a, a, and 8. For simplicity we will work in 
the alignment limit 8 = a + 7/2, where the free param- 
eters are reduced to 6, ay, with tan 0, = — cot 8 tan ap, 
and the 125 GeV Higgs boson is SM-like, if there is no 
CP-violation. The overall contribution to the eEDM is 
then 
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with f"(z,%) = —8(af(z) + g(z)/x)/3) and g'(z,2) = 
(3f(z) + 5g(z))x. The contributions from neutral Higgs 
with the Z gauge boson and charged Higgs with W gauge 
boson as the propagator are generally smaller, so we ne- 
glect them in the calculation [I9]. 
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yy ww* LL" 
ATLAS|1.17 + 0.27 220]] 0.997935 1.441035 
CMS | 1.144038 0.727070 0.937032 
bb TT 
ATLAS] 0.52 + 0.40 iA 
CMS [1.15 Æ 0.62 0.78 + 0.27 


TABLE I: The LHC data used for the fitting. 
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FIG. 2: The allowed parameter region |œp| versus tan 8 from 
the Higgs global fits and EDM experiments. The grey, yellow, 
green region are the 1,2,3 o allowed region for the LHC Higgs 
fits. The left and right region bounded by the black dashed 
lines and the upper region bounded by the green solid line are 
excluded by the eEDM and neutron EDM experiment at 90% 
C.L. The constraint from Mercury experiment is weaker 
than the one from neutron EDM so we do not plot it. The 
light blue region is the theoretical bound from the Higgs mass 
spectra Mp = 125 GeV, m p+ =420 GeV, my2 = 400 GeV and 
my3 = 450 GeV Bl. The blue solid contours represent the 
CP-phase of the top Yukawa coupling tan 6;. The orange band 
and contours give the baryon-to-entropy density ratio in unit 
of 107". To calculate the BAU, a wall velocity vw = 0.02, 
a wall width Lù = 5/T. and a critical temperature Te = 
100 GeV are assumed for the bubbles generated during the 
EWPT. 


The fitting results are presented in Fig. [2] where the 
inclusive LHC data published in March 2013 (see Ta- 
ble |I) and the most recent ACME results [12] are ap- 
plied. In the presence of CP-even and CP-odd Higgs mix- 
ing, both the current ACME constraints and the Higgs 
global fits favors the region with tan 8 ~ 1, where we have 
tan, ~ — tan ap. The former is easy to understand since 
a cancellation between f’(z;,tan 8) and g’(zw,tan 8) in 
Eq. (15) requires tan@ = 1.04. The latter is because 
a relatively small tan 6 can help avoid a too large signal 
rate of h — bb, and hence an over-suppressed h > yy rate 
(see Table [I). In this cancellation region, a CP-violation 
effect with |tanap|,|tan6@,| > 0.1 is allowed while the 
most stringent constraints are from the nEDM. 


MSSM 


Though the MSSM is of type II 2HDM, there is no 
tree-level CP-violation in the Higgs sector either explic- 
itly or spontaneously [58], due to a vanishing A; term 
in Eq. (7). So, the CP-phases used for EWBG mostly 


arise in the tree-level superparticle sectors, such as the 
chargino, neutrolino, squark and slepton sectors. 


The explicit CP-violation in these sectors can break 
the CP symmetry in the Higgs sector at loop level, lead- 
ing to CP-even and CP-odd mixing terms in the Higgs 
squared mass matrix. But the Higgs CP-mixture caused 
by this effect is small due to the loop suppression. For 
nonstandard Higgs bosons we notice that the CP-mixture 
is typically below 10%, consistent with [32], even if the 
CP-phase arises from the stop sector, while for the SM- 
like Higgs boson, the CP-mixture is suppressed more by 
an extra tan 6 factor. So the Higgs eigenstate are approx- 
imately CP-eigenstates, with their couplings with elec- 
trons are either |c.| >> |če| (for CP-even Higgs bosons) 
or |ce| < |é| (for CP-odd Higgs bosons). On the other 
hand, a relatively large tan 8 is favored in the MSSM, 
given that the tree-level mass of the SM-like Higgs boson 
is larger in this case. This leads to || < |cH4], and 
hence a small h contribution to the eEDM (mainly via 
the c.é, term). So in the MSSM with the EWBG mecha- 
nism implemented, the main contributions to the eEDM 
are made by nonstandard Higgs bosons unless they are 
highly decoupled. 


Among these CP-violating sources, the one arising in 
the chargino sector is of particular interest because of its 
high efficiency in generating the BAU via EWBG. The 
charginos enter the Hyy and the Ayy loops as new medi- 
ators, inducing non-trivial contributions to the eEDM via 
the ceé, and ¢.c, terms, respectively. Though these two 
contributions are comparable in magnitude, their signs 
are typically different. Given an extra minus sign for the 
term cĉ, in Eq. (5), this scenario is strongly constrained 
by the ACME eEDM bound [6]. However, recall that 
charged particles like stau leptons enter the Ayy loop as 
well. If a non-trivial CP-phase is turned on in the stau 
sector, new contributions to the eEDM will be introduced 
via the čec, term. (Note, such a CP-violating coupling 
will not induce non-trivial contributions to the eEDM 
via the Hyy loop or the ceé, term, since stau leptons 
are scalar particles.) This provides a potential cancel- 
lation, such that a CP-phase in the Higgs-chargino cou- 
plings which is large enough for implementing the EWBG 
mechanism, is still allowed. 


With the CP-violation turned on in the stau sector, 
the main contributions to the eEDM are then given by 
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with the terms in the first and the second lines mediated 


by A and H, respectively. Here C = ##5, z2” = T2 and 
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As for mg+, €+(CP-odd interaction between H and 


charginos), we refer to [33] for the details. 
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FIG. 3: Parameter region allowed by the current eEDM, 
mercury EDM and neutron EDM bounds, and favored by the 
EWBG (orange band) in the u — tan 8 plane. The chargino, 
stau contributions to the eEDM are plotted with black dashed 
and blue dashed contours, respectively. The magenta and red 
regions have been excluded by the ACME and the Mercury 
experiments |30|. The bounds of neutron EDM are much 
weaker and fall outside of the figure. Here we are assum- 
ing: charged Higgs mass my,+ = 450 GeV; trilinear softly 
SUSY-breaking parameters (except A+) 1.5 TeV and Ay = 
2.5 TeV; soft masses of gauginos Mı = 0.1 TeV, M2 = 0.45 
TeV and M3 = 3.5 TeV; soft masses of squarks and slep- 
tons 0.1MQ,,0,,p, = 0.1Ma@z2,t2,D2 = Ma3,u3,p3 = 1.2 TeV, 
and 0.1Mz, = 0.1Mr, = Mr, = 0.15 TeV, 0.05Mp, = 
0.05Mz, = Me, = 0.35 TeV, and CP-phases Arg(uMž) = 
Arg(uA;) = 90°. With their mass set to be ~ 10 TeV, the 
contributions of the squarks and sleptons of the first two gen- 
erations to the one-loop EDM are suppressed. To calculate 
the baryon asymmetry generated via EWBG, we closely follow 
the work done in [31], with an expansion velocity vw = 0.02, 
a wall width L = 5/Te and a critical temperature Te = 100 
GeV assumed for the bubbles generated during the EWPT. 


Fig. [3] depicts all low energy experimental constraints 
calculated by the CPsuperH [33] [34]. Note, however, the 
CPsuperH codes used in the analysis was revised by the 
authors or this letter (for details, see Appendix C), in 
which we corrected a sign error in the anomalous dimen- 
sion of the dipole operators, updated the QCD hadron 
matrix elements, and incorporated the mixing effects be- 
tween the dipole and psedo-scalar operators during renor- 
malization group running, and the missed contributions 


to the eEDM of the W-mediated hyy loop in the Barr- 
Zee diagrams. With these corrections, we notice that the 
neutron and mercury bounds may become substantially 
weaker than the ones given by the original CPsuperH. 

As indicated in Fig. the charginos have a neg- 
ative contribution to the eEDM (black-dashed con- 
tours), by coupling with both A and H. The staus, 
on the other hand, have a positive contribution to 
the eEDM (blue-dashed contours), by mainly coupling 
with A. Both contributions are enhanced by tan 8 
because of G4 œ tang and cH œ 1/cos6. Their 
dependences on the u parameter however is differ- 
ent. os charginos, the H/A-chragino- chee cou- 
plings a ro i(Cr)i1(Cr) 39/2 — h.c. and Tae ge 
Opal? — h.c., i = 1,2 [B5]. With u > Mə as- 
sumed here, a small u sale wills increase the off-diagonal 
term in the chargino mixing matrices, (C',)12 and (CR)ai, 
and hence the overall eEDM contribution. As for staus, 
a larger u leads to a lighter 7, because of a larger mixture 
term x |wtan$—A,|. Therefore, their contribution to 
the eEDM increases for a larger u value. Due to their 
cancellation there exists a blank region in Fig. [3] where 
the total eEDM is below the current ACME bound. This 
region overlaps with the EWBG favored region which has 
been excluded by the ACME bound, if only the chargino 
contribution is taken into account. One benchmark point 
is presented in Table 


Ryy tan 8 p(TeV) mi, (GeV) 
0.84 12.8 1.31 125.4 
[de/e] (cm) [dirg/e](cm) |dn/el(em) m/s 


4.2x10 7 3.0x10 7 9.6x10~7" 0.85x10~™ 


TABLE II: An benchmark in the MSSM, with the other pa- 
rameters set as F igb] 


DISCUSSION AND CONCLUSION 


In this letter, we present a class of cancellation mech- 
anisms to suppress the total contributions of Barr-Zee 
diagrams to the eEDM, which may occur either between 
the contributions of a CP-mixing Higgs boson, with the 
other Higgs bosons being decoupled, or between the con- 
tributions of a CP-even and a CP-odd Higgs bosons. 
As an illustration, we study two scenarios: the type- 
I 2HDM where the tree-level CP-phase arises from the 
Higgs sector, and the MSSM where the tree-level CP- 
phase arises from Higgs-superparticle interaction sectors. 
In the 2HDM, tan ~ 1 is favored by the LHC Higgs 
bounds and the contributions of the Barr-Zee diagrams 
to the eEDM are mainly mediated by a CP-mixed SM- 
like Higgs boson. With a cancellation between them, a 
CP-phase as large as O(0.1 — 1) is still allowed for the 


top Yukawa coupling, induced by the Higgs CP mixing. 
In the MSSM, a large tan 8 is favored by the LHC Higgs 
bounds, and the contributions of the Barr-Zee diagrams 
to the eEDM are mainly mediated by nonstandard CP- 
even and CP-odd neutral Higgs bosons (if they are not 
as heavy as 10 TeV scale or above). With a cancellation 
between them, a maximal CP-phase in their couplings 
with superparticles like charginos is still allowed. In both 
cases a successful EWBG is possible, which turns out to 
be challenging without the assistance of these cancella- 
tion mechanisms. 

The EWPT, another key element for the EWBG, is 
not addressed in this letter, since it is “orthogonal” to 
the discussions above. It is straightforward to generalize 
the discussions on the eEDM cancellation to the 2HDM 
+ a singlet, or the MSSM extensions with an extra singlet 
superfield or gauge group, where a strong enough EWPT 
is not difficult to achieve (e.g., see [B6H38]. Furthermore, 
in the MSSM while an extremely light right handed stop 
scenario (~ 110 GeV) has been ruled out by the 
Higgs global fits [41}/43] and direct searches, a moder- 
ately light stop with another light scalar particles (like 
sbotoom or stau) may still be viable [44]. The detailed 
discussions in this regard will be left for future studies. 
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Neutron and Mercury EDMs in the MSSM: 
the CPsuperH vs. Our Analysis 


In the Appendix, we will introduce the modifications 
and corrections made to the CPsuperH in our analysis. 
Then we will make a comparison between the constraints 
of neutron and Mercury EDMs given by the original 
CPsuperH and the modified one. As a start, let’s simply 
introduce the theoretical methods used for calculating 
the neutron and Mercury EDMs, basically following 


ref BA Bd]. 


(1) Renormalization Group Running of Wilson 
Coefficients 


The EDM of a fermion (d7; electron EDM is an excep- 
tion which is denoted by de in this letter), the Chromo- 


EDM of a quark (dF), and the Weinberg operator cou- 
pling (d°) are defined by the Lagrangian 
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Here F” and G*#” are the electromagnetic and strong 
field strengths, T* = \*/2 denotes the generators of the 
SU(3)c group. In the MSSM, the Weinberg operator dS 
is given by: 


d@ = (a° Han (19) 


with (d@)9 being contributed by the quark-squark-gluino 
coupling [45] and (dC)! being contributed by neutral 


Higgs bosons [46] [47]. 


Defining the wilson coefficients as 
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with my and A denoting the fermion masses and the 
MSSM scale respectively, the Lagrangian can be 
rewritten as 
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To calculate the neutron and mercury EDMs, we need 
to incorporate the effect of renormalization group run- 
ning of the WIlson coefficients from new physics scale to 
a hadron scale. During this process, flavor-conserving 
CP-odd four-fermion operators may lead to nontrivial 
corrections to the Wilson coefficients of the CEDM and 
Weinberg operators via mixing. So a more complete La- 
grangian for the calculation of the neutron and mercury 
EDMs should be 


Lopv = L(c)epM + 3 - oi 


Here the first two CP-odd four-fermion operators 
oi 
oF 


qqq 1754, 
q'd'qiysq (23) 


can be generated through CP-violating neutral Higgs- 
boson mixing in the t-channel and CP-violating Yukawa 
threshold corrections. The corresponding CP-odd coeffi- 
cients are given by 


Cg Č 
ci = Ia Ia WIR’ 
xg Cg' Č 
CT? = Iq' Iq M ’ (24) 


with gq(q’) = Mqcq’)/V and v = 2Mw /g. The last CP-odd 
four-fermion operator 


OF Oe Go” dG iO pv Y5la » (25) 


on the other hand, is generated from the operatror mixing 
effects of C44 and C4 which follow the Eq. (27) below. 


To calculate the Wilson coefficients (= ô -39 ) 


A a) 
at a GeV scale, we need to take m evolution 
for C = (# jp 292,04, 644, Ga", Cf) from the 
MSSM scale A down to GeV scale, based on the Renor- 
malization Group Equations (RGE) : 


d 
din uu 


C=C T (26) 


The one-loop anomalous dimension matrix is given by 


Eys (0) 0 


i ie Uy Ysf IYF 0 , (27) 
ay? yp 0 ae 
with 
+8CPr 0 0 
Ys = |+8Cr +16Cr — 4N 0 3 (28) 
0 +2N N +2nf + Bo 
vp = [-12Cr +6], (29) 
—12Cpr 0 —1 
yp=| 0 -12Cr —1 , (30) 
12 12 -8Cr- $ 
0 0 0 
Yaf = , Q’ ae 0 , (31) 
1672 5% 167e 0 
and 
yep = [+4 +4 0]. (32) 
Here N = 3, Cr = (N? —1)/(2N) = 4/3, bo = (11N — 


2n,y)/3, ng is the flavor number, and q runs over u, d, b,. 


Besides the RG running and operator mixing effects 
calculated from Eq. (26) which give corrections to the 
Wilson coefficients ôg, ôg and Ce, there are two addi- 
tional bottom quark threshold effects which needs to be 
emphasized below. First, the Weinberg operator receives 
a shift from the bottom quark CEDM after the bottom 
quark is integrated out at m, [BI 53]: 


(33) 


Here 44(ms) = ô? (mp) + Adz (ms) is the b-quark CEDM 
at the mp scale, with 5°(mp) being the direct CEDM 
calculated from two-loop Barr-Zee graphs. The Ady (mp) 
is the bottom quark CEDM correction from RG running 
and operator mixing. Keeping the leading logarithmic 
terms that contributes to Aĝp(m,) at the matching scale 
L = Mp, we have 

Aae S aCO$(Mn) igg u o (89 
which could be figured out from Eq. and is from 
integrating out b quark at one-loop level. 


The second important corrections by integrating out 
the bottom quark are the shift to CEDMs of quarks: 


we g Mp 
Nia 2 
6474 mg 
: ‘ Mn \* 
(C(t) + ČP tm) (10g 2) (385) 
b 


which are actually a two loop effects where the operator 
Cy! and CY’ first mix into Č and then mix into 
Adq(mp). Numerically, the suppression from one addi- 
tional loop is compensated by the m,/m, enhancement. 


(2) Neutron and Mercury EDMs 
The neutron EDM is calculated by 
dn = (eC#5y + C454) + (eið J ecdõa) + BECe (36) 


In the CPsuperH, we updated hadronic matrix elements 
with C = 0.82 x 1078, C= —3.3 x 1078, ce = 0.82 x 
1078, ¢4 = 1.63 x 1078 and 67 = 2 x 107% ecm 54. 

Compared to the electron EDM de and the CP-odd 
electron-nucleon interactions 


L= CsēiyseN N 


+ Cpēe Niys N + Ohee Niys N, (37) 


the nuclear Schiff moment (S) has a larger contribution 
to the mercury EDM. The Schiff moment is generated by 
long-range, pion-exchange mediated P- and T-violating 
nucleon-nucleon interactions, 


PINY =Ñ [arr t+ gi) n° + gP (237° —T-7)| N38) 


In a general context, the isoscalar and isovector cou- 
plings gt, gr’ 
gP [54], so the mercury EDM is approximately given 


by [54], 


are dominant over the isotensor coupling 


2 
dug = kg S & rs" (a03 T age) , (39) 
in which 
GL = ño (bu + ba) + 1%)Ca: 
IP = fay (bu — ba) + 1G)Ce | 
ga © 1.26, F, = 186 MeV. (40) 


To do the calculation, in the CPsuperH we updated 
hadronic matrix elements with [54]: ij) = —2 x 1077, 


fia) = —4 x 1077, 1) ~ VG) = 2x 107°; updated 


the nuclear matrix elements with [55]: ao = 0.01 efm, 
ay = +0.02 e fm’; and assume a new atomic sensitivity 


coefficient xg = —2.8 x 1074 fm~? [55]. At last we would 
like to emphasize again that both the sign correction of 
the anomalous dimension coefficient ye and the mixing 
effect of the RGE operators have been incorporated in 
calculating {ôu, ôa, Ôu, Ôd, CG} at a hadron scale in our 
analysis, which may cause an important change for the 
bounds of both the neutron and mercury EDMs given 
by the original CPsuperH. 


(3) The CPsuperH vs. Our Analysis 


In the CPsuperH, we made the following corrections 
and modifications: 


1. A sign error in an anomalous-dimension coefficient 
Ye was fixed, which is equal to the (1,1) compo- 
nent of Eq. (28), up to some numerical factor 
caused by different definitions. The CPsuperH fol- 
lows Ref. where ye = 8/3 is defined. As pointed 
in Ref. [48] [50], however, ye should be —8/3 rather 
than 8/3 in the notation of Ref. [49]. 


2. The mixing effects between the Weinberg’s gluonic 
operator, quark color dipole operator, quark elec- 
tric dipole operator and the CP-odd four-fermion 
operators which are defined by Eq. are in- 
corporated in the CPsuperH. 


3. The hadronic matrix elements are updated. 


4. The Barr-Zee diagrams with the Higgs-photon- 
photon(Z boson) loop mediated by a W boson are 
incorporated. 


These corrections and modifications made to the 
CPsuperH may significantly change the neutron and the 
mercury EDM bounds for CP-violation in the MSSM. As 
a comparison, the neutron and mercury EDM bounds af- 
ter and before the revision of the CPsuperH are shown in 


Fig. 


and Fig. |4| respectively, with a chiral model is 


assumed for hadrons and the same parameter values as- 
sumed in both plots. The comparison indicates that both 
the neutron and mercury EDM bounds are overestimated 
by the original CPsuperH for the assumed parameter val- 


ues. 
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FIG. 4: The neutron and mercury EDM bounds before the 
revision of the CPsuperH, with the parameter values assumed 
to be the same as the ones defined in Fig.[3] The green and the 
red regions are excluded by the neutron and mercury EDM 
bounds, respectively. 
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